Maternal circadian rhythm disruption affects neonatal inflammation via metabolic reprograming of myeloid cells.

Disruption of circadian rhythm during pregnancy produces adverse health outcomes in offspring; however, the role of maternal circadian rhythms in the immune system of infants and their susceptibility to inflammation remains poorly understood. Here we show that disruption of circadian rhythms in pregnant mice profoundly aggravates the severity of neonatal inflammatory disorders in both male and female offspring, such as necrotizing enterocolitis and sepsis. The diminished maternal production of docosahexaenoic acid (DHA) and the impaired immunosuppressive function of neonatal myeloid-derived suppressor cells (MDSCs) contribute to this phenomenon. Mechanistically, DHA enhances the immunosuppressive function of MDSCs via PPARγ-mediated mitochondrial oxidative phosphorylation. Transfer of MDSCs or perinatal supplementation of DHA relieves neonatal inflammation induced by maternal rhythm disruption. These observations collectively demonstrate a previously unrecognized role of maternal circadian rhythms in the control of neonatal inflammation via metabolic reprograming of myeloid cells.

Evaluation of the efficacy and safety of immunotherapy in sarcoma: a two-center study.

Background: Sarcoma is a highly heterogeneous malignancy with a poor prognosis. Although chemotherapy and targeted therapy have improved the prognosis to some extent, the efficacy remains unsatisfactory in some patients. The efficacy and safety of immunotherapy in sarcoma need further evaluation. Methods: We conducted a two-center study of sarcoma patients receiving PD-1 immunotherapy at Tianjin Medical University Cancer Institute and Hospital and Henan Provincial Cancer Hospital. The treatment regimens included PD-1 inhibitor monotherapy and combination therapy based on PD-1 inhibitors. The observed primary endpoints were median progression-free survival (mPFS) and median overall survival (mOS). Survival curves were compared using the Kaplan-Meier method. Results: A total of 43 patients were included from the two centers. The median follow-up time for all patients was 13 months (range, 1-48 months). In the group of 37 patients with advanced or unresectable sarcoma, the mPFS was 6 months (95%CI: 5-12 months), and the mOS was 16 months (95%CI: 10-28 months). The ORR was 10.8% (4/37), and the DCR was 18.9% (7/37). Subgroup analysis showed no significant differences in mPFS (p=0.11) and mOS (p=0.88) between patients with PD-L1 negative/positive expression. There were also no significant differences in mPFS (p=0.13) or mOS (p=0.72) between PD-1 inhibitor monotherapy and combination therapy. Additionally, there were no significant differences in mPFS (p=0.52) or mOS (p=0.49) between osteogenic sarcoma and soft tissue sarcoma. Furthermore, the results showed no significant differences in mPFS (p=0.66) or mOS (p=0.96) between PD-1 inhibitors combined with targeted therapy and PD-1 inhibitors combined with AI chemotherapy. Among the 6 patients receiving adjuvant therapy after surgery, the mPFS was 15 months (95%CI: 6-NA months), and the mOS was not reached. In terms of safety, most adverse events were mild (grade 1-2) and manageable. The most severe grade 4 adverse events were bone marrow suppression, which occurred in 4 patients but resolved after treatment. There was also one case of a grade 4 adverse event related to hypertension. Conclusion: Immunotherapy is an effective treatment modality for sarcoma with manageable safety. Further inclusion of more patients or prospective clinical trials is needed to validate these findings.

SingleQ: a comprehensive database of single-cell expression quantitative trait loci (sc-eQTLs) cross human tissues.

Mapping of expression quantitative trait loci (eQTLs) and other molecular QTLs can help characterize the modes of action of disease-associated genetic variants. However, current eQTL databases present data from bulk RNA-seq approaches, which cannot shed light on the cell type- and environment-specific regulation of disease-associated genetic variants. Here, we introduce our Single-cell eQTL Interactive Database which collects single-cell eQTL (sc-eQTL) datasets and provides online visualization of sc-eQTLs across different cell types in a user-friendly manner. Although sc-eQTL mapping is still in its early stage, our database curates the most comprehensive summary statistics of sc-eQTLs published to date. sc-eQTL studies have revolutionized our understanding of gene regulation in specific cellular contexts, and we anticipate that our database will further accelerate the research of functional genomics. Database URL: http://www.sqraolab.com/scqtl.

Target genes regulated by CLEC16A intronic region associated with common variable immunodeficiency.

BACKGROUND: CLEC16A intron 19 has been identified as a candidate locus for common variable immunodeficiency (CVID). OBJECTIVES: This study sought to elucidate the molecular mechanism by which variants at the CLEC16A intronic locus may contribute to the pathogenesis of CVID. METHODS: The investigators performed fine-mapping of the CLEC16A locus in a CVID cohort, then deleted the candidate functional SNP in T-cell lines by the CRISPR-Cas9 technique and conducted RNA-sequencing to identify target gene(s). The interactions between the CLEC16A locus and its target genes were identified using circular chromosome conformation capture. The transcription factor complexes mediating the chromatin interactions were determined by proteomic approach. The molecular pathways regulated by the CLEC16A locus were examined by RNA-sequencing and reverse phase protein array. RESULTS: This study showed that the CLEC16A locus is an enhancer regulating expression of multiple target genes including a distant gene ATF7IP2 through chromatin interactions. Distinct transcription factor complexes mediate the chromatin interactions in an allele-specific manner. Disruption of the CLEC16A locus affects the AKT signaling pathway, as well as the molecular response of CD4+ T cells to immune stimulation. CONCLUSIONS: Through multiomics and targeted experimental approaches, this study elucidated the underlying target genes and signaling pathways involved in the genetic association of CLEC16A with CVID, and highlighted plausible molecular targets for developing novel therapeutics.

An interpretable model predicts visual outcomes of no light perception eyes after open globe injury.

BACKGROUND: The visual outcome of open globe injury (OGI)-no light perception (NLP) eyes is unpredictable traditionally. This study aimed to develop a model to predict the visual outcomes of vitrectomy surgery in OGI-NLP eyes using a machine learning algorithm and to provide an interpretable system for the prediction results. METHODS: Clinical data of 459 OGI-NLP eyes were retrospectively collected from 19 medical centres across China to establish a training data set for developing a model, called 'VisionGo', which can predict the visual outcome of the patients involved and compare with the Ocular Trauma Score (OTS). Another 72 cases were retrospectively collected and used for human-machine comparison, and an additional 27 cases were prospectively collected for real-world validation of the model. The SHapley Additive exPlanations method was applied to analyse feature contribution to the model. An online platform was built for real-world application. RESULTS: The area under the receiver operating characteristic curve (AUC) of VisionGo was 0.75 and 0.90 in previtrectomy and intravitrectomy application scenarios, which was much higher than the OTS (AUC=0.49). VisionGo showed better performance than ophthalmologists in both previtrectomy and intravitrectomy application scenarios (AUC=0.73 vs 0.57 and 0.87 vs 0.64). In real-world validation, VisionGo achieved an AUC of 0.60 and 0.91 in previtrectomy and intravitrectomy application scenarios. Feature contribution analysis indicated that wound length-related indicators, vitreous status and retina-related indicators contributed highly to visual outcomes. CONCLUSIONS: VisionGo has achieved an accurate and reliable prediction in visual outcome after vitrectomy for OGI-NLP eyes.

Maternal antibiotic exposure enhances ILC2 activation in neonates via downregulation of IFN1 signaling.

Microbiota have an important function in shaping and priming neonatal immunity, although the cellular and molecular mechanisms underlying these effects remain obscure. Here we report that prenatal antibiotic exposure causes significant elevation of group 2 innate lymphoid cells (ILC2s) in neonatal lungs, in both cell numbers and functionality. Downregulation of type 1 interferon signaling in ILC2s due to diminished production of microbiota-derived butyrate represents the underlying mechanism. Mice lacking butyrate receptor GPR41 (Gpr41-/-) or type 1 interferon receptor IFNAR1 (Ifnar1-/-) recapitulate the phenotype of neonatal ILC2s upon maternal antibiotic exposure. Furthermore, prenatal antibiotic exposure induces epigenetic changes in ILC2s and has a long-lasting deteriorative effect on allergic airway inflammation in adult offspring. Prenatal supplementation of butyrate ameliorates airway inflammation in adult mice born to antibiotic-exposed dams. These observations demonstrate an essential role for the microbiota in the control of type 2 innate immunity at the neonatal stage, which suggests a therapeutic window for treating asthma in early life.

Unraveling the causal genes and transcriptomic determinants of human telomere length.

Telomere length (TL) shortening is a pivotal indicator of biological aging and is associated with many human diseases. The genetic determinates of human TL have been widely investigated, however, most existing studies were conducted based on adult tissues which are heavily influenced by lifetime exposure. Based on the analyses of terminal restriction fragment (TRF) length of telomere, individual genotypes, and gene expressions on 166 healthy placental tissues, we systematically interrogate TL-modulated genes and their potential functions. We discover that the TL in the placenta is comparatively longer than in other adult tissues, but exhibiting an intra-tissue homogeneity. Trans-ancestral TL genome-wide association studies (GWASs) on 644,553 individuals identify 20 newly discovered genetic associations and provide increased polygenic determination of human TL. Next, we integrate the powerful TL GWAS with placental expression quantitative trait locus (eQTL) mapping to prioritize 23 likely causal genes, among which 4 are functionally validated, including MMUT, RRM1, KIAA1429, and YWHAZ. Finally, modeling transcriptomic signatures and TRF-based TL improve the prediction performance of human TL. This study deepens our understanding of causal genes and transcriptomic determinants of human TL, promoting the mechanistic research on fine-grained TL regulation.

Landscape of enhancer disruption and functional screen in melanoma cells.

BACKGROUND: The high mutation rate throughout the entire melanoma genome presents a major challenge in stratifying true driver events from the background mutations. Numerous recurrent non-coding alterations, such as those in enhancers, can shape tumor evolution, thereby emphasizing the importance in systematically deciphering enhancer disruptions in melanoma. RESULTS: Here, we leveraged 297 melanoma whole-genome sequencing samples to prioritize highly recurrent regions. By performing a genome-scale CRISPR interference (CRISPRi) screen on highly recurrent region-associated enhancers in melanoma cells, we identified 66 significant hits which could have tumor-suppressive roles. These functional enhancers show unique mutational patterns independent of classical significantly mutated genes in melanoma. Target gene analysis for the essential enhancers reveal many known and hidden mechanisms underlying melanoma growth. Utilizing extensive functional validation experiments, we demonstrate that a super enhancer element could modulate melanoma cell proliferation by targeting MEF2A, and another distal enhancer is able to sustain PTEN tumor-suppressive potential via long-range interactions. CONCLUSIONS: Our study establishes a catalogue of crucial enhancers and their target genes in melanoma growth and progression, and illuminates the identification of novel mechanisms of dysregulation for melanoma driver genes and new therapeutic targeting strategies.

ADAR1 links R-loop homeostasis to ATR activation in replication stress response.

Unscheduled R-loops are a major source of replication stress and DNA damage. R-loop-induced replication defects are sensed and suppressed by ATR kinase, whereas it is not known whether R-loop itself is actively involved in ATR activation and, if so, how this is achieved. Here, we report that the nuclear form of RNA-editing enzyme ADAR1 promotes ATR activation and resolves genome-wide R-loops, a process that requires its double-stranded RNA-binding domains. Mechanistically, ADAR1 interacts with TOPBP1 and facilitates its loading on perturbed replication forks by enhancing the association of TOPBP1 with RAD9 of the 9-1-1 complex. When replication is inhibited, DNA-RNA hybrid competes with TOPBP1 for ADAR1 binding to promote the translocation of ADAR1 from damaged fork to accumulate at R-loop region. There, ADAR1 recruits RNA helicases DHX9 and DDX21 to unwind R-loops, simultaneously allowing TOPBP1 to stimulate ATR more efficiently. Collectively, we propose that the tempo-spatially regulated assembly of ADAR1-nucleated protein complexes link R-loop clearance and ATR activation, while R-loops crosstalk with blocked replication forks by transposing ADAR1 to finetune ATR activity and safeguard the genome.

Inferring CTCF-binding patterns and anchored loops across human tissues and cell types.

CCCTC-binding factor (CTCF) is a transcription regulator with a complex role in gene regulation. The recognition and effects of CTCF on DNA sequences, chromosome barriers, and enhancer blocking are not well understood. Existing computational tools struggle to assess the regulatory potential of CTCF-binding sites and their impact on chromatin loop formation. Here we have developed a deep-learning model, DeepAnchor, to accurately characterize CTCF binding using high-resolution genomic/epigenomic features. This has revealed distinct chromatin and sequence patterns for CTCF-mediated insulation and looping. An optimized implementation of a previous loop model based on DeepAnchor score excels in predicting CTCF-anchored loops. We have established a compendium of CTCF-anchored loops across 52 human tissue/cell types, and this suggests that genomic disruption of these loops could be a general mechanism of disease pathogenesis. These computational models and resources can help investigate how CTCF-mediated cis-regulatory elements shape context-specific gene regulation in cell development and disease progression.

Cross-ancestry genome-wide association meta-analyses of hippocampal and subfield volumes.

The hippocampus is critical for memory and cognition and neuropsychiatric disorders, and its subfields differ in architecture and function. Genome-wide association studies on hippocampal and subfield volumes are mainly conducted in European populations; however, other ancestral populations are under-represented. Here we conduct cross-ancestry genome-wide association meta-analyses in 65,791 individuals for hippocampal volume and 38,977 for subfield volumes, including 7,009 individuals of East Asian ancestry. We identify 339 variant-trait associations at P < 1.13 × 10-9 for 44 hippocampal traits, including 23 new associations. Common genetic variants have similar effects on hippocampal traits across ancestries, although ancestry-specific associations exist. Cross-ancestry analysis improves the fine-mapping precision and the prediction performance of polygenic scores in under-represented populations. These genetic variants are enriched for Wnt signaling and neuron differentiation and affect cognition, emotion and neuropsychiatric disorders. These findings may provide insight into the genetic architectures of hippocampal and subfield volumes.

Noncoding SNP at rs1663689 represses ADGRG6 via interchromosomal interaction and reduces lung cancer progression.

A previous genome-wide association study (GWAS) revealed an association of the noncoding SNP rs1663689 with susceptibility to lung cancer in the Chinese population. However, the underlying mechanism is unknown. In this study, using allele-specific 4C-seq in heterozygous lung cancer cells combined with epigenetic information from CRISPR/Cas9-edited cell lines, we show that the rs1663689 C/C variant represses the expression of ADGRG6, a gene located on a separate chromosome, through an interchromosomal interaction of the rs1663689 bearing region with the ADGRG6 promoter. This reduces downstream cAMP-PKA signaling and subsequently tumor growth both in vitro and in xenograft models. Using patient-derived organoids, we show that rs1663689 T/T-but not C/C-bearing lung tumors are sensitive to the PKA inhibitor H89, potentially informing therapeutic strategies. Our study identifies a genetic variant-mediated interchromosomal interaction underlying ADGRG6 regulation and suggests that targeting the cAMP-PKA signaling pathway may be beneficial in lung cancer patients bearing the homozygous risk genotype at rs1663689.

Systematic fine-mapping and functional studies of prostate cancer risk variants.

To date, genome-wide association studies (GWAS) have revealed over 200 genetic risk loci associated with prostate cancer; yet, true disease-causing variants remain elusive. Identification of causal variants and their targets from association signals is complicated by high linkage disequilibrium and limited availability of functional genomics data for specific tissue/cell types. Here, we integrated statistical fine-mapping and functional annotation from prostate-specific epigenomic profiles, 3D genome features, and quantitative trait loci data to distinguish causal variants from associations and identify target genes. Our fine-mapping analysis yielded 3,395 likely causal variants, and multiscale functional annotation linked them to 487 target genes. We prioritized rs10486567 as a genome-wide top-ranked SNP and predicted HOTTIP as its target. Deletion of the rs10486567-associated enhancer in prostate cancer cells decreased their capacity for invasive migration. HOTTIP overexpression in enhancer-KO cell lines rescued defective invasive migration. Furthermore, we found that rs10486567 regulates HOTTIP through allele-specific long-range chromatin interaction.

GBC: a parallel toolkit based on highly addressable byte-encoding blocks for extremely large-scale genotypes of species.

Whole -genome sequencing projects of millions of subjects contain enormous genotypes, entailing a huge memory burden and time for computation. Here, we present GBC, a toolkit for rapidly compressing large-scale genotypes into highly addressable byte-encoding blocks under an optimized parallel framework. We demonstrate that GBC is up to 1000 times faster than state-of-the-art methods to access and manage compressed large-scale genotypes while maintaining a competitive compression ratio. We also showed that conventional analysis would be substantially sped up if built on GBC to access genotypes of a large population. GBC's data structure and algorithms are valuable for accelerating large-scale genomic research.

Genome-wide analysis of genetic pleiotropy and causal genes across three age-related ocular disorders.

Age-related macular degeneration (AMD), cataract, and glaucoma are leading causes of blindness worldwide. Previous genome-wide association studies (GWASs) have revealed a variety of susceptible loci associated with age-related ocular disorders, yet the genetic pleiotropy and causal genes across these diseases remain poorly understood. By leveraging large-scale genetic and observational data from ocular disease GWASs and UK Biobank (UKBB), we found significant pairwise genetic correlations and consistent epidemiological associations among these ocular disorders. Cross-disease meta-analysis uncovered seven pleiotropic loci, three of which were replicated in an additional cohort. Integration of variants in pleiotropic loci and multiple single-cell omics data identified that Müller cells and astrocytes were likely trait-related cell types underlying ocular comorbidity. In addition, we comprehensively integrated eye-specific gene expression quantitative loci (eQTLs), epigenomic profiling, and 3D genome data to prioritize causal pleiotropic genes. We found that pleiotropic genes were essential in nerve development and eye pigmentation, and targetable by aflibercept and pilocarpine for the treatment of AMD and glaucoma. These findings will not only facilitate the mechanistic research of ocular comorbidities but also benefit the therapeutic optimization of age-related ocular diseases.

PR-DUB safeguards Polycomb repression through H2AK119ub1 restriction.

Polycomb group (PcG) proteins are critical chromatin regulators for cell fate control. The mono-ubiquitylation on histone H2AK119 (H2AK119ub1) is one of the well-recognized mechanisms for Polycomb repressive complex 1 (PRC1)-mediated transcription repression. Unexpectedly, the specific H2AK119 deubiquitylation complex composed by additional sex comb-like proteins and BAP1 has also been genetically characterized as Polycomb repressive deubiquitnase (PR-DUB) for unclear reasons. However, it remains a mystery whether and how PR-DUB deficiency affects chromatin states and cell fates through impaired PcG silencing. Here through a careful epigenomic analysis, we demonstrate that a bulk of H2AK119ub1 is diffusely distributed away from promoter regions and their enrichment is positively correlated with PRC1 occupancy. Upon deletion of Asxl2 in mouse embryonic stem cells (ESCs), a pervasive gain of H2AK119ub1 is coincident with increased PRC1 sampling at chromatin. Accordingly, PRC1 is significantly lost from a subset of highly occupied promoters, leading to impaired silencing of associated genes before and after lineage differentiation of Asxl2-null ESCs. Therefore, our study highlights the importance of genome-wide H2AK119ub1 restriction by PR-DUB in safeguarding robust PRC1 deposition and its roles in developmental regulation.

An autoimmune pleiotropic SNP modulates IRF5 alternative promoter usage through ZBTB3-mediated chromatin looping.

Genetic sharing is extensively observed for autoimmune diseases, but the causal variants and their underlying molecular mechanisms remain largely unknown. Through systematic investigation of autoimmune disease pleiotropic loci, we found most of these shared genetic effects are transmitted from regulatory code. We used an evidence-based strategy to functionally prioritize causal pleiotropic variants and identify their target genes. A top-ranked pleiotropic variant, rs4728142, yielded many lines of evidence as being causal. Mechanistically, the rs4728142-containing region interacts with the IRF5 alternative promoter in an allele-specific manner and orchestrates its upstream enhancer to regulate IRF5 alternative promoter usage through chromatin looping. A putative structural regulator, ZBTB3, mediates the allele-specific loop to promote IRF5-short transcript expression at the rs4728142 risk allele, resulting in IRF5 overactivation and M1 macrophage polarization. Together, our findings establish a causal mechanism between the regulatory variant and fine-scale molecular phenotype underlying the dysfunction of pleiotropic genes in human autoimmunity.

eaQTLdb: An atlas of enhancer activity quantitative trait loci across cancer types.

Enhancers are key regulatory elements that exert crucial roles in diverse biological processes, including tumorigenesis and cancer development. Active enhancers could produce transcripts termed enhancer RNAs (eRNAs), which could be used as an index of enhancer activity. Here, we present a versatile data portal, enhancer activity quantitative trait loci database (eaQTLdb; http://www.bioailab.com:3838/eaQTLdb), for exploring the effects of genetic variants on enhancer activity and prioritizing candidate variants across different cancer types. By leveraging the accumulated multiomics data, we systematically identified genetic variants which influence enhancer activity in different cancer types, termed as eaQTLs. We have linked the eaQTLs to hallmarks of cancer and patients' overall survival to illustrate their potential biological roles in cancer development and progression. Notably, eaQTLs associated with the infiltration abundance of 24 different immune cell types were identified and incorporated into eaQTLdb. In addition, we applied colocalization analyses to examine 59 complex diseases and traits to identify eaQTLs colocalized with diseases/traits GWAS signals. Overall, eaQTLdb, incorporating a rich resource for exploration of eaQTLs in different cancer types, will not only benefit users in prioritizing candidate genetic variants and enhancers, but also help researchers decipher the roles of eaQTLs in the dysregulated pathways of cancer and tumor immune microenvironment, opening new diagnostic and therapeutic avenues in precise medicine.

Retrograde regulation of mitochondrial fission and epithelial to mesenchymal transition in hepatocellular carcinoma by GCN5L1.

Metabolic reprogram is crucial to support cancer cell growth and movement as well as determine cell fate. Mitochondrial protein acetylation regulates mitochondrial metabolism, which is relevant to cancer cell migration and invasion. The functional role of mitochondrial protein acetylation on cancer cell migration remains unclear. General control of amino acid synthesis 5 like-1(GCN5L1), as the regulator of mitochondrial protein acetylation, functions on metabolic reprogramming in mouse livers. In this study, we find that GCN5L1 expression is significantly decreased in metastatic HCC tissues. Loss of GCN5L1 promotes reactive oxygen species (ROS) generation through enhanced fatty acid oxidation (FAO), followed by activation of cellular ERK and DRP1 to promote mitochondrial fission and epithelia to mesenchymal transition (EMT) to boost cell migration. Moreover, palmitate and carnitine-stimulated FAO promotes mitochondrial fission and EMT gene expression to activate HCC cell migration. On the other hand, increased cellular acetyl-CoA level, the product of FAO, enhances HCC cell migration. Taken together, our finding uncovers the metastasis suppressor role as well as the underlying mechanism of GCN5L1 in HCC and also provides evidence of FAO retrograde control of HCC metastasis.

CD73, a Promising Therapeutic Target of Diclofenac, Promotes Metastasis of Pancreatic Cancer through a Nucleotidase Independent Mechanism.

CD73, a cell surface-bound nucleotidase, facilitates extracellular adenosine formation by hydrolyzing 5'-AMP to adenosine. Several studies have shown that CD73 plays an essential role in immune escape, cell proliferation and tumor angiogenesis, making it an attractive target for cancer therapies. However, there are limited clinical benefits associated with the mainstream enzymatic inhibitors of CD73, suggesting that the mechanism underlying the role of CD73 in tumor progression is more complex than anticipated, and further investigation is necessary. In this study, CD73 is found to overexpress in the cytoplasm of pancreatic ductal adenocarcinoma (PDAC) cells and promotes metastasis in a nucleotidase-independent manner, which cannot be restrained by the CD73 monoclonal antibodies or small-molecule enzymatic inhibitors. Furthermore, CD73 promotes the metastasis of PDAC by binding to the E3 ligase TRIM21, competing with the Snail for its binding site. Additionally, a CD73 transcriptional inhibitor, diclofenac, a non-steroidal anti-inflammatory drug, is more effective than the CD73 blocking antibody for the treatment of PDAC metastasis. Diclofenac also enhances the therapeutic efficacy of gemcitabine in the spontaneous KPC (LSL-KrasG12D/+ , LSL-Trp53R172H/+ , and Pdx-1-Cre) pancreatic cancer model. Therefore, diclofenac may be an effective anti-CD73 therapy, when used alone or in combination with gemcitabine-based chemotherapy regimen, for metastatic PDAC.